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The response of mesoscopic superconductors to an ac magnetic field is inves- 
tigated both experimentally and with numerical simulations. We study small 
square samples with dimensions of the order of the penetration depth. We 
obtain the ac susceptibitity X = x' + ^x" c-^ microwave frequencies as a func- 
tion of the dc magnetic field Hdc- We find that the dissipation, given by x" > 
has a non monotonous behavior in mesoscopic samples. In the numerical 
simulations we obtain that the dissipation increases before the penetration of 
vortices and then it decreases abruptly after vortices have entered the sample. 
This is verified experimentally, where we find that x" has strong oscillations 
as a function of Hd^. in small squares of Pb. 
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The response of superconductors to an ac magnetic field has been of 
interest in the last year^. The microwave surface impedance and the ac 
magnetic susceptibility X = x' + ^x" have been extensively studied. The 
imaginary part of the susceptibility, x" : is proportional to the dissipation in 
the sample. In macroscopic type II superconductors, the ac dissipation in- 
creases with magnetic field, proportional to the vortex density. Very recently, 
there has been an interest in the study of mesoscopic superconductos where 
the sample dimensions are of the order of the London penetration deptlP. In 
this work we will show that in mesoscopic superconductors the ac dissipation 
has a non monotonous and oscillating dependence with magnetic field. 

We first show the results of numerical simulations of mesoscopic squares 
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Fig. 1. ac dissipation, 'x" , and number of vortices, Ny, vs. magnetic field H. 
(a) Large sample, 40A x 40A. (b) Mesoscopic sample, lOA x lOA 



using the time-dependent Ginzburg-Landau (TDGL) equationsP^ 

'^[{V - iA)H + [l - T){1 - \^^)^ (1) 

r)Im[^'*(V -iA)^'] - K^V X V X A (2) 

where ^ and A are the order parameter and vector potential respectively 
and T is the temperature. Lengths are normalized by the coherence length 
^(0), times by to = 47rcr„Ai/c^ A by i7c2(0)^(0), T by Tc; and we have 
T] = /{Airani^'^D). The boundary conditions are (V — iA)-^"^ = and 
Bz\s = (V X A)z\s = H':/^. We consider an ac+dc magnetic field: Hf^* = 
Hdc + hac cos{ujt) with hac ^ Hdc- The magnetization is given by 47rM(t) = 
{B[t)) — H{t), with {B(t)) averaged over the sample. The ac magnetic 
susceptibilities are obtained as x' = J^^ M{t) cos{u;t)d{ujt) and x" = 
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Jq^ M{t) sm{Lut)d{ujt). We consider here k = 2, T = 0.5 and the pa- 
rameter T] = 12. We use a finite difference discretization schem^^ with 
Ax = Ay = 0.5^ and time step At < O.OlSio- Fig. 1(a) shows x'' {Hdc) 
at u; = 0.06fo obtained in a "large" sample of 40A x 40A. We also show 
in Fig. 1(a) the number of vortices as a function of Hdc- At small Hdc, in 
the Meissner state, x" increases continuously with increasing Hdc- In this 
case, the magnetic field induces static supercurrents which deplete the or- 
der parameter at the boundary and this leads to an increase of dissipation 
for increasing H^c- Above a penetration field Hp > Hd (there is a sur- 
face barrieJ^ vortices enter into the sample and the dissipation increases 
linearly with the number of vortices. In Fig. 1(b) we show the similar case 
but for a mesoscopic square of size lOA x lOA. When the magnetic field is 
increased above the first penetration field Hp = Hpi, the first vortices enter 
into the sample, and we see that there is a discontinuous jump in x" with 
a decrease of dissipation just at Hp. Further jumps in x" i^H^f^ are present 
at the other magnetic fields for vortex penetration, Hp2, Hp^, . . ., where x" 
decreases. The jumps are followed by a later continuous increase of x"{H(ic) 
with increasing H^c while the number of vortices remains fixed. This 
increase of dissipation at fixed has been interpreted as due to "nascent 
vortices"'^^ which give a large fluctuation of |^'| at the boundary. Later, when 
vortices enter, there is a sudden decrease of the ac losses. 

We now discuss the experimental measurements. An array of 200 x 250 
Pb squares were deposited by thermal evaporation on a silicon subtrate. The 
squares of L = 1.3 ^m on side and d = 300A thick where defined by e-bem 
lithography and a subsequent lift-off. This thickness d corresponds to an 
effective Kg// = Kp^^Pb/d ~ 1.3 and to a length L ?a lOAg//. A SEM image 
of part of the sample is shown in the lower inset of Fig. 2. Magnetic res- 
onance measurements have been made with a commercial Bruker ESP 300 
spectrometer. Spectra were acquired at 9.5 GHz (X-Band) using a TE102 
rectangular cavity. The external magnetic field was modulated with a 100 
kHz ac field of varying amplitude to obtain the field derivative of the mi- 
crowave absorption {dx" / dH). Low temperature measurements were made 
using a continuous fiow helium cryostat. The sample was positioned inside 
the cavity with the dc field and the rf field perpendicular and parallel to the 
film plane, respectively. Spectra features were generally of very weak inten- 
sity and it was necessary to use long acquisition times (several minutes) in 
order to obtain a reasonably good signal to noise ratio. In this way, spectra 
of dx" /dH vs. H were obtained. Fig. 2 shows the integrated curve, which is 
proportional to x" ^ function of the the magnetic field H for a temperature 
T = A.bK. We first observe an overall increase of x" with H which is also 
observed in bulk Ph and it is due to microscopic mechanisms (not described 
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Fig. 2. Measured ac susceptibility, Y cx x"; in 1-3 Pb squares at 9.5 GHz. 
Lower inset: SEM image of part of the sample. 

by the TDGL equations). An enhancement of the behavior of x" at low 
fields is shown in the upper inset of Fig. 2. Wc sec that the ac dissipation 
has a non monotonous oscillatory behavior as a function of H. This is in 
qualitative agreement with the theoretical result of Fig. 1(b). Effects of dis- 
order due to small variations in the shape of the squares may cause that the 
oscillations of x" arc not as uniform as in Fig. 1(b). In conclusion, we have 
observed that in mesoscopic squares the dependence of ac dissipation with 
magnetic field is very sensitive to vortex entrance events. 
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